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In continuation of [19], we discuss the question of spectral mini- 

p ^ mal 3-partitions for the rectangle ] — |, |[x] — |, |[, with < a < b. 
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3-partition is obtained by the three nodal domains of the third eigen- 
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mechanism of transition for increasing ^ between these nodal mini- 

j> mal 3-partitions and non nodal minimal 3-partitions at the value \/ 1 

and discuss the existence of symmetric candidates for giving minimal 

00 3-partitions when < f < 1- Numerical analysis leads very natu- 

rally to nice questions of isospectrality which are solved by introducing 
Aharonov-Bohm Hamiltonians or by going on the double covering of 

QQ the punctured rectangle. 
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analysis. In the same spirit, we discuss here the similar question for the 
rectangle Ila,b •=]~t;i[^]~|;|[; 'with < a < 6. For a given partitioij^ D 
of an open set Q hj k open subsets Di, we can consider 

A{V) = max A(A) , (1-1) 

i=l,...,k 

where X{Di) is the ground state energy of the Dirichlet Laplacian on Di . We 
denote the infimum on every ^-partitions of Q by 

= inf A(V) . (1.2) 
We look for minimal /^-partitions, that is partitions such that iifc(^) = A{V) . 



It has been observed in [IS] that, when < | < y | the minimal 3- 
partition is given by the three nodal domains of the third eigenfunction cor- 

2'2l- -l 2' el-'J 2'2^ J 6'6l 



responding to the three rectangles ] — ^, ^[x] — ^, — ] — ^, ^[x] — ^, M and 



a al^lb b\ 
2 ' 2 I- J 6 ' 2 I 



In the case when < f < 1 , we can show as for the disk, see [9], that 
Sl3{Ra,b) is not an eigenvalue. Indeed, A2(Ra,6) = A3(Ra,6) < A4(Ra,fe)- By 



Theorem 2.4 below, see [19] for the proof, £3(Ra^b) cannot be an eigenvalue 



and hence the associated minimal partition cannot be nodal. 



We will describe in Section |4] a possible mechanism of transition for in- 
creasing I between these nodal minimal 3-partitions and non nodal minimal 

3-partitions at the value ^J~^ and discuss the existence of symmetric candi- 
dates for giving minimal 3-partitions when < f < 1 in Sections 6 7 

We can exhibit numerically some candidates for the minimal 3-partition 
using symmetry. Assuming that there is a minimal partition which is sym- 
metric with respect to the axis {y = 0} , and intersecting the partition with 
the half-square ] — ^,^[x]0,|[, one is reduced to analyze a family of Dirichlet- 
Neumann problems. Numerical computation^ performed by V. Bonnaillie- 
Noel and G. Vial (in January 2006) lead to a natural candidate V for a 
symmetric minimal partition (see Figure [T(a)[ ). We observe numerically that 
the three lines of N{T>) (i.e. the interior boundary of the subdomains Di of 



the partition, see Definition (2.6)) meet at the center (0, 0) of the square. As 



^See the next section for precise definitions. 

^See http;//w3. bretagne.ens-cachan.fr/math/simulations/MinimalPartitions/ 
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(a) First candidate T) (b) Second candidate I?""^™ 

Figure 1: Candidates for the minimal 3-partition of the square. 



expected by the theory they meet at this critical point with equal angle ^ 
and meet the boundary orthogonally. This choice of symmetry is not unique. 
By exploring numerically the possibility of a minimal partition with another 
symmetry (diagonal), we get the surprise to find another candidate T)^*^^ 
with A(P"^'") = A{V) (see Figure p^(b)| . 



This leads very naturally to nice questions of isospectrality which are solved 
using Aharonov-Bohm Hamiltonians or by going on the double covering of 
the punctured rectangle. Sections [5]j6] concern these questions. 



The paper is organized as follows: 
In Section [2| we recall some notation and properties concerning the nodal 
partitions. 

In Section [3} we start with the analysis of the rectangle based on results of 
[19] and enumerate in particular all the possible Courant-sharp situations. 
There appears the limiting case a/b = y^3/8 detailed in Section [i] in which 
a mechanism is proposed which explains the transition between the nodal 
3-partition and the non nodal one. 

Motivated by the numerical simulations for the square we find two candidates 
for the minimal 3-partition with the same energy. We analyze in Section [5] 
the Aharonov-Bohm Hamiltonian and give some isospectral properties for 
rectangles. This theory is applied in Section |6] to explain numerical simula- 
tions for the minimal 3-partitions on the square. The paper ends in Section [7] 
with some heuristics on the deformation of symmetric minimal partitions cor- 
roborated with numerical simulations for rectangles ] — 



2 ' 2 



TV TT I 
2 ' 2 I 



from 



to e = 1. 
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2 Definitions, notation and previous results 



As in fL9] (see also [20]), we consider the Dirichlet Laplacian on a bounded 
domain C M^, which is piecewise C°°. We denote, for any open domain D, 
the lowest eigenvalue of the Dirichlet realization H{D) of —A in D by A(D). 
For any function u G Cq{Q) , we introduce 



N{u) = {xen\ u{x) = 0} (2.1) 

and call the components of i7 \ N{u) the nodal domains of u. The number 
of nodal domains of such a function is denoted by ^{u). 

We now recall the main definitions and results concerning spectral min- 
imal partitions and refer to [19] for proofs and details. For k > 1 {k & N), 
we call a fc-partition of Q a family V = {-Di}f=i of pairwise disjoint open 
domains in Q. Such a partition is called strong if 



Int (Uti A) \dn = n. (2.2) 



We denote by Ok the set of such partitions. 
For V G Ok we introduce 



A{V) = max A(A) , (2.3) 

i=l,...,k 



and 



£fc = inf A(V) . (2.4) 

feOfe 

V is called a (spectral)]^ minimal fc-partition if £a;(^) = A(D) and a nodal 
minimal fc-partition if V consists of the nodal domains of an eigenfunction 

To each strong partition V we associate a graph G{T>) in the following way: 
We say Z)i, Dj G D are neighbors, and we denote this by Di ~ Dj, if 



Int (A U Dj) \ dVt is connected. (2.5) 

We associate to each Di & V a. vertex Vi and for each pair Di ~ Dj an 
edge Cij. This defines a planar graph G{T>). We say that the partition is 
admissible if the corresponding graph is bipartite. We recall that a nodal 
partition is always admissible. 

Attached to a partition we can associate a closed set G defined by 



N{V) = [j{dDi n n) . (2.6) 

i 

^We will omit the word spectral. 
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This leads us to introduce the set Ai{Q) of regular closed sets which share 
with nodal sets all the standard properties except at isolated critical points 
where they have onljj^the equal angle meeting property. More precisely, 
we recall from [19] the following definition. 

Definition 2.1 

A dosed set N G Q belongs to Ai{fl) if N satisfies 

(i) There are finitely many distinct Xi G Q (1 N and associated positive 
integers z/(a;j) > 3 such that, in a sufficiently small neighborhood of 
each Xi, N is the union of ulxi) smooth arcs (non self- crossing) with 
one end at Xi and such that in the complement of these points in Vt, 
N is locally diffeomorphic to a regular curve. The set of these critical 
points of N is denoted by X{N). 

(a) dQ n N consists of a finite set of points Zi such that at each Zi, p{zi) 
arcs hit the boundary with p{zi) > 1. We denote the set of critical 
points of Nndn by Y{N) . 

A partition V is called regular if the corresponding N{T>) is regular. Let us 
now recall the main theorems. 

Theorem 2.2 

For any k there exists a minimal regular strong k-partition and any minimal 
k-partition admits a representative which is regular and strong. 

The existence of a minimal regular strong partition has been shown|^ by 
Conti-Terracini- Verzini in [TTl [T31 [12] , while the second part of the theorem 
has been shown in jl9] . 

In the following, we always consider the regular representative without men- 
tioning it explicitly. We have now the following converse theorem, see [T9] . 

Theorem 2.3 

Assume that there is a minimal admissible k-partition. Then this partition 
is associated to the nodal set of an eigenfunction corresponding to 2k{^). 

'^We do not assume anymore that the number of hues arriving at a critical point is 
even. 

^See also [ig and HQj. 
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This result was completed in [19j in relation with the Courant-sharp property. 
We recall that if u is an eigenfunction of the Dirichlet Laplacian in Q attached 
to the k-th eigenvalue A^, then Courant's Theorem says that the number of 
nodal domains fi{u) satisfies ^{u) < k . Pleijel's Theorem says that, when 
the dimension is > 2, then the previous inequality is strict for k large. 
As in [2], we say that u is Courant-sharp if fi{u) = k. For any integer 
/c > 1, we denote by Lk{Q) the smallest eigenvalue for which the eigenspace 
contains an eigenfunction with k nodal domains. In general we have 

Xk{n) < < Lfc(fi) . (2.7) 

The next result of [19] gives the full picture of the equality cases: 

Theorem 2.4 

//£fc((]) = Lfc(fi) or Xk{n) = S.k{^), then 

In addition, any minimal (regular) k-partition is a nodal partition corre- 
sponding to an eigenfunction associated to Xk{^)- 

As a consequence of Euler's Formula, we have described in [20] the possi- 
ble topological types of a non admissible minimal 3-partition of a connected 
regular open set. 

Proposition 2.5 

Let Q be simply connected and let us consider a minimal 3-partition T) = 
{Di, D2, D-i) of Q associated to £3(^2) . Let us suppose that 

XsiQ) < . (2.8) 

For N = N{V), we denote by vi^Xi) and p{zi) the number of arcs associated 
with Xi e X{N), respectively Zi G Y{N^. Then there are three possibilities 
(see Figure^: 

(a) X{N) consists of one point x with z/(x) = 3 and Y{N) consists of either 
three distinct yi, 1/2, ys points with p{yi) = p{y2) = pivs) = 1 , two distinct 
points yi, 1/2 with p{yi) = 2, p{y2) = I or one point y with p{y) = 3 . 

(b) X{N) consists of two distinct points Xi, X2 with z/(xi) = i'{x2) = 3 . 
Y{N) consists either of two points yi, y2 such that p{yi) + p(l/2) = 2 or of 
one point y with p{y) = 2 . 

(c) X{N) consists again of two distinct points xi, X2 with ^{xi) = i'{x2) = 3, 
but Y{N) = . 
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Figure 2: The three configurations (a), (b) and (c), with an additional sym- 
metry with respect to the x-axis. 



3 A first analysis of the rectangle 

This part is taken from [19]. Note that when is a rectangle Ra^h := 
] ~ t'i[x] ~ |;|[; spectrum of H(Iia,b) and the properties of the eigen- 
functions are analyzed as toy models in |26i §4]. The spectrum is given by 



Am,n := ttM ^ + ^ ) with (m, n) G {Wf . 



These eigenvalues are simple if is irrational. Except for specific remarks 
concerning the square or the case when a? /IP' = 3/8, we assume 

is irrational. (3.1) 

So we can associate to each eigenvalue Xm,n, up to a non-zero multiplicative 
constant, a unique eigenfunction Um,n such that fi{um,n) = fnn. Given k G 
N*, the lowest eigenvalue corresponding to k nodal domains is given, at least 



under Assumption (3.1), by 



2 2 

m n 



Lfc(R,,,) = vr^ inf ( ^ + ^ ) . (3.2) 

mn=k y 0^ / 

In the case when is rational we could have problems in the case of 

multiplicities. We have then to analyze a continuous families of nodal sets 
of eigenfunctions living in an eigenspace of dimension > 1. We will see in 
Section |4] that it is just for these values that new nodal partitions may appear, 
which could be, by deformation, the starting point of non admissible minimal 
partitions. 

We now recall all the possible Courant-sharp situation^ described in 



^We do not know whether for certain rational additional Courant sharp eigen- 

values could show up. 
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, 3 5 

(i) m = 3, n = 2 and - < -n^ < - ■ 

3 

(ii) m = 2, n = 2 and - < 777 < 1 . 

5 0^ 

. , 3 
(m) m = 1, n > 1 and — < 



62 - n2 - 1 



If we now focus on the case A; = 3, we get that A3(Ra^b) is Courant-sharp iff 
a^/fe^ < 3/8. Hence, the hmiting situation is 

a2 _ 3 
62 ~ 8 ■ 

This corresponds to a double eigenvalue and to the pairs (m, n) = (1, 3) and 
(m,n) = (2,l). 



4 Transition from Courant-sharp to a non nodal 
minimal partition 

We start from a rectangle with a = vre and b = tt and would like to analyze 
£3(6) := £3(R^e,7r). The critical situation corresponds to 

6=^378- (4.1) 

So the first result (deduced from [19]) which was recalled in the previous 
section writes: 

Proposition 4.1 

(z) Ife < ^/3/8, then 2.3(e) = 9 + l/e^ and £3(6) is an eigenvalue. 
(11) If v/378 < e < 1, then £3(6) < 9 + l/e^ . 

We would like to understand how the transition occurs at e = ^378. We 
make the assumption that in the deformation the minimal partition remains 
symmetric with respect to y = 0. This is indeed the case for e < ^3/8 , 
because the eigenfunction corresponding to Ai^3 is cos | cos 3?/ and the cor- 
responding nodal lines are composed of two horizontal lines y = —71/6 and 
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y = 7r/6. This is also the case for e = 

because all the eigenf unctions 
have this symmetry and any minimal partition is nodal. 



The numerical computations for the square (see Figure 1(a) and [9", §3.2]) 
push us to conjecture that the nodal lines N{T>) for the minimal 3-partition 
V is for e > -^/S/S is the union of a segment [— vr/2, xo(e)] (on the line |/ = 0) 
and of two symmetric arcs connecting the point (xo(e),0) to the boundary 
of the rectangle (up and down). 



The first conjecture is that xo(e) is increasing monotonically from — 7r/2 
to for e g] a/3/8, 1]. This has been partially verified (admitting the sym- 
metry and that the minimal 3-partition is of type (a)) numerically. 
The second conjecture is that the minimal 3-partition will "tend" as e tends 
to a/3/8 from above to a nodal partition, losing there its non bipartite char- 
acter. 

The point here is that, when e = a/3/8, we have an eigenvalue of multi- 
plicity two giving the possibility of constructing a continuous family of nodal 
minimal 3-partitions. For this, we consider the family 

X 2x 
'^a,i3{x, y) = a cos — cos 3y + /3 sin — cosy , 

with + 7^ 0, and analyze their zero set. Of course, for t 7^ 0, ipa,ii and 
ipta,tf} have the same zero set. 



We first show that the zero set of ipa,p has no critical point inside the 
rectangle for the critical value of e. Using the factorization of (pa,i3 in the 
form 

(Pa,/3{x, y) = COS y cos - ya [l — 4 sin^ y) + 2(3 sin — j , 
we observe that ipa,i3 = is equivalent inside the rectangle to ipa,i3 = with 

ipa^p{x,y) := a(l - isin^y) + 2/3 sin - . (4.2) 



Hence we now look at the (closure of the) zero set of tpa,(3 in the rectangle 
and particularly to the critical points inside and at the boundary. If we look 
at the zeros of dxilJa,i3, we obtain 



X 

j3 cos - = . 
e 
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This implies /5 = and we get that (pa,o{x, y) = a cos 3y cos ^ . Hence we get 
that, for any {a, (3) 7^ (0,0), there is no critical point inside this rectangle. 
It remains to look at what is going on at the boundary and to determine the 
singular points where two lines touch. An analysis of the function {x, y) 1— > 
a(l — 4sin^?/) + 2/? sin - at the boundary shows that critical points at the 
boundary can only occur for y = and a ± 2/5 = 0. 

Hence we have obtained that the only nodal sets having critical sets are (up 
to a multiplicative constant) the nodal domains of the eigenf unctions ip2,i 
and f2-i- 

Figure [3] gives the nodal set of the functions (fa,/3 for several values of (a,/3). 




5 Aharonov-Bohm Hamiltonian and isospec- 
trality 

As explained in the introduction, this new analysis is motivated by numerical 
computations showing that pushing the same idea as in [9J but using the 
symmetry with respect to the diagonal, one gets the same eigenvalue and 
again a 3-partition with unique singular point at the center. This will be 



further explained in more detail in Subsection 6.1 
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5.1 Basic material 



This material appears already in JX7] and is motivated by the work of Berger- 
Rubinstein If f2 is an open set such that e a possibility is to 
consider the Aharonov-Bohm Laplacian in the punctured Q = Q \ {0}, with 
the singularity of the potential at the center and normalized flux $ = 1/2. 
The magnetic potential with flux $ takes the form 

A{x,y) = iA,ix,y),A2ix,y)) = ^ ("4' 4) • (^-l) 

We know that the magnetic field vanishes identically in and, in any cut 
domain (such that it becomes simply connected), one has 

Axdx ^ A2dy = ^ dQ , (5.2) 

where 

z = X + iy = r expiO . (5.3) 

So the Aharonov-Bohm operator in any open set 17 C \ {0} will always 
be defined by considering the Friedrichs extension starting from C^{Cl) and 
the associated differential operator is 

- Aa := {D, - Ai)2 + {Dy - A^f . (5.4) 

In polar coordinates (which of course are not very well adapted to the square 
but permit a good analysis at the origin), the Aharonov-Bohm Laplacian 
reads: 

,^ 1 sin^^,, , 1 cos 
- Aa = {D^ + - —f + {Dy - - —f , (5.5) 

or 

d"^ I d I 1 2 

This operator is preserving "real" functions in some modified sense. Follow- 
ing [T7], we will say that a function u is K-real, if it satisfies 

Ku = u, (5.7) 

where K is an antilinear operator in the form 

K = expi9T, (5.8) 

and where F is the complex conjugation 

Tu = u. (5.9) 
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The fact that — Aa preserves K-real eigenfunctions is an immediate conse- 
quence of 

Ko{-Aj^) = i-AA)oK. (5.10) 



Remark 5.1 

Note that our choice of K is not unique. = expiaK is also antilinear, 



satisfies (5.10) and 



Ki = ld. 



As observed in [T7j, it is easy to find a basis of K-real eigenfunctions. These 
eigenfunctions (which can be identified with real antisymmetric eigenfunc- 
tions of the Laplacian on a suitable double covering VL^ of VL) have a nice 
nodal structure (which is locally in the covering the same as the nodal set 
of real eigenfunctions of the Laplacian), with the specific property that the 
number of lines in Vt ending at the origin should be odd. More generally a 
path of index one around the origin should always meet an odd number of 
nodal lines. 



5.2 Symmetries of the rectangle 

We consider now a domain Vl which has the symmetries of a rectangle. More 
precisely, if we denote by ai and (T2 the symmetries respectively defined by 

(^i{.x,y) = {-x,y) , 0-2(0;, y) = {x,-y) , (5.11) 

we assume that 

ain = h , a2Cl = n. (5.12) 
We assume that Q is convex (to simplify) and write 

nn{y = o} =]-^,|[x{0}, 

and 

nn{x = o} = {o}x]-^-,^-[. 

If El is the natural action on L'^{Cl) associated with cxi 

^Mx,y) = u{-x,y) , (5.13) 
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we observe that the Aharonov-Bohm operator does not commute with Si 
but with the antihnear operator 

S5:=irSi. (5.14) 

So if u is an eigenfunction, E^m is an eigenfunction. 

Moreover, and this explains the choice of "z" before FSi, since K and 

commute, 

Koj:i = j:iok, (5.15) 

S^M is also a K-real eigenfunction if m is a K-real eigenfunction. One can do 
the same thing with E2, associated with (T2, 

T.2u{x,y) = u{x,-y) . (5.16) 

This leads this time to 

S^ = rS2. (5.17) 

Similarly, we have 

/s:oS^ = S^oi^ , (5.18) 

hence if u is a K-real function, is also a K-real eigenfunction. 
We now show 



Proposition 5.2 

If Cl has the symmetries of the rectangle (5.11), then the multiplicity of the 
groundstate energy of —Aa is 2. 

More generally the multiplicity of any eigenvalue is even. 



Proof: 

As observed in , we can reduce the analysis of the Aharonov-Bohm Hamil- 
tonian to the ii'-real space where 

Lj^{tl) = {ue L'^itl) , Ku = u}. 

The scalar product on Lj^, making of Lj^ a real Hilbert space, is obtained by 
restricting the scalar product on L'^{Cl) to L|- and it is immediate to verify 
that {u , v) is indeed real for u and v in 
Observing now that 

S^oS^ = /, (5.19) 

we obtain by writing 

u = \{I + T.l)u+\{I-^l)u, 
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an orthogonal decomposition of Lj^ into 

^x = ^x,s.©^Ls,, (5.20) 

where 

^x,si = {ueL^K ^ ^lu = u }, 

and 

We have just to show that the restriction Hi of ^(/ + S^) to Lj^ 

1 
2 



ni:=-(/ + S5)/^|^, (5.21) 



is a projector. It is indeed clear that Hi is (]R-)linear and that 11^ = Hi . It 
remains to verify that 11^ = Hi. But we have, for u, v in Lj^, 

(S^M , v) = i{Tv , Sim) = ^(rSif , u) = (S^t> , u) = {u , S^t>) . 



Moreover the decomposition (5.20) is respected by — Aa. 



Similarly, one can define the projector 112 by restriction of S2 to Lj^. 



The second statement of Proposition 5.2 will be a consequence of the follow- 
ing lemma 

Lemma 5.3 

Let 

= , (5.22) 
then S3 commutes with — Aa and := (Sg)/^!^ is a unitary operator from 
^Isi onto 

Proof: We note that 

E=3 = zE3, (5.23) 

where S3 is associated with: 

a3{x,y) = i-x,-y). (5.24) 

The lemma follows then from the property that if m is a solution of Ku = u 
and S^M = M, then 

S^ S^ M = S^ S^ S^ M = -S^ S^ S^ M = -S^ u , 

where we have used the anticommutation of S^ and Sg: 

SJ S^ = -S^ SJ . (5.25) 

□ 
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It remains to show that the first eigenvalue has multiphcity 2. We already 
know that it has an even multiplicity. It is enough to prove that the mul- 
tiplicity is at most 2. Here we can use the results of |17]. Actually those 
results have to be extended slightljj^ since they are obtained assuming that 
the domain is homeomorphic to an annulus (so we are in a limiting case). It 
has been shown in [T7] that the nodal set of a K-real groundstate is a line 
joining the center to the outer boundary. If the multiplicity of the ground- 
state eigenvalue is strictly greater than 2 we can, as in [T7j, construct by 
linear combination of eigenfunctions a groundstate for which the zeroset hits 
the outer boundary at two distinct points, hence a contradiction. □ 

We observe that the proof of the proposition gives more explicitly the 
decomposition of — Aa on into the direct orthogonal sum of two unitary 
equivalent Hamiltonians. What we have done with can be similarly done 
with Eg. This gives immediately the following proposition. 

Proposition 5.4 

The four following operators -Aa,Si, -AA,aSi, -Aa,E2 and -AA,aS2 re- 
spectively defined by the restriction of — Aa to L\y.^, L\^y,^, L^^.^^^ and 
Lj^aT,2 are isospectral to — Aa- Moreover A is an eigenvalue of any of the 
first four operators with multiplicity k{X) if and only if X is an eigenvalue of 
multiplicity 2k{X) of — Aa ■ 

Now we would like to analyze the nodal patterns of eigenfunctions in the 
various symmetry spaces. 

Lemma 5.5 

Ifu G C°°{Cl) (^Lj^s^ ^hen its nodal set contains [— 1,0] x {0}. Moreover, in 
fi\{]-f,0]x{0}}, 

id 

V = exp 'a 

satisfies 

Tv = V and II2V = v . 

Proof: 

Noting that Ku = u and = u for y = and x < we immediately 
'^The paper of Alziary-Fleckinger-Takac [T] is dealing with this case. 
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obtain that u{x, 0) = u{x, 0) = —u{x, 0). Hence u{x, 0) = for x < 0. 
An immediate computation gives 

ie ^ ie ie 

Vv — exp —Vu — exp — — Ku — exp —— u — v , 



and 



2^21' = exp — L2M = exp — Vu — V . 



□ 



Wc would like to compare to some Dirichlet-Neumann problems on half- 
domains. We call upperhalf-r2 the set 



Q"'^ = Q n {y > 0} , 



(5.26) 



and introduce similarly the lowerhalf, lefthalf and righthalf domains defined 

by 

n^'* = n < 0} , i^'"^ = n < 0} , f)"'' = n {x > 0} . (5.27) 

The previous lemma leads to 



Dirichlot 

_i 



Neumann 





<- Dirichlet 


Qleh 






<- Neumann 



Figure 4: Domains Q"'', Q"', Q'^'^, Q"'*. 



Proposition 5.6 

If u is a K-real E2 invariant eigenfunction of — Aa; then the restriction to 
Quh ^jg-jjp_^^ Q j^gQ^ eigenfunction of the realization of the Laplacian 
in Q^'^, with the following Dirichlet-Neumann condition at dfl^^: Dirichlet 
except on ]0, |] x {0} where we put Neumann. 

In particular, if X is an eigenvalue of —A a, then A is an eigenvalue of the 
Laplacian in D,^^ with this Dirichlet-Neumann condition. 
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Conversely, if we consider an eigenfunction u of the Laplacian in with 
this Dirichlet-Neumann condition and extend it into m*^^* by symmetry in 
nW - f,0] X {0}, then 

V = u^""' (5.28) 
is a K-real eigenfunction of the Aharonov-Bohm Laplacian. More precisely, 



the function v is first defined by the formula (5.28) for 6' g] — vr, tt[ and then 
extended as a function on the punctured square. Due to the properties 
of u, one can verify that v is in the form domain of the Aharonov-Bohm 
operator. Starting from a groundstate of the DA^-problem in Q^^, we get a 
ii'-real eigenstate in Lj^ of the Aharonov-Bohm operator. 

One can do the same type of argument in the three other half-domains 
and obtain the following result. 

Proposition 5.7 

The following problems have the same eigenvalues: 

• The Dirichlet problem for the Aharonov-Bohm operator on Vl. 

• The Dirichlet-Neumann problem for the Laplacian on VT^ . 

• The Dirichlet-Neumann problem for the Laplacian on Q''^^. 

• The Dirichlet-Neumann problem for the Laplacian on fi'^ . 

• The Dirichlet-Neumann problem for the Laplacian on 

Of course this applies in particular to the case of the rectangle. 

Let us go a little further by giving explicitly a unitary operator from 
-^Icsi onto -Z^lcsa proving the isospectrality. This is the object of 



Lemma 5.8 

The operator 



t^2i:=^(/ + S^) 



is a unitary operator from L\ onto L\ j,^ whose inverse is given by 

^12:=^(/ + ^). 
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Proof: Let u G then, using (5.25), 

UMiu = ^(/ + + + ElE',)u = ^(/ + + - E^EDu = u . 

The proof that U21U12 = I on L\-y,2 obtained in the same way. Let us 
prove the norm conservation. If m G L\ , then 



\U2^u 



But if Tj\u = u, we can write 



This leads to 



\U21U 



\u\\\yu G 



(5.29) 
□ 



6 Application to minimal 3-partitions 
6.1 Discussion on the square 

We look at the first excited eigenvalue of the Dirichlet problem in the punc- 
tured square. The rules of [T7] give constraints about the nodal structure of 
the ii'-real eigenfunctions, which were already used. In particular we should 
have an odd number of lines arriving at the center. So it is clear that {0} 
belongs to the nodal set. If three lines arrive at and if the nodal partition 
is a 3-partition of type (a) it is reasonable to hope that this will give us a 
minimal 3-partition. 



Let us explain the numerical strategy developed in [9', §3] to exhibit a 
candidate for the minimal 3-partition of the square. According to Theo- 
rem 2.3, if the minimal 3-partition of the square is admissible, it is associatd 



to the nodal set of an eigenf unction for A3. But there is no such function 
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and therefore the minimal 3-partition is non bipartite. Then we look for non 



bipartite 3-partitions whose topologies are enumerated in Proposition 2.5 



and illustrated by Figure [2j We first use the axial symmetry along the axis 
{y = 0}. To recover a partition of type (a), (b) or (c), we compute the second 
eigenfunction and the next ones of the mixed Dirichlet-Neumann Laplacian 
in with Dirichlet conditions except respectively on 



• [xo,a/2] X {0} for type (a); 

• [xo,Xi] X {0} for type (b); 

• [-a/2,xo] X {0} U [xi,a/2] x {0} for type (c). 



These boundary conditions are illustrated in Figure [5j We move the points xq 
and Xi along the segment [—a/2, a/2] x {0}. We expect to find an eigenfunc- 
tion such that, after symmetrization, their associated nodal sets constitute 
a 3-partition and the nodal lines meet at the interior critical point with an 
angle of 27r/3. 





Figure 5: Mixed problems with Dirichlet-Neumann conditions. 



We recall the guess inherited from numerical computations. When mini- 
mizing over Dirichlet-Neumann problems in Q'-^ (by putting Dirichlet except 
Neumann on ]xq, |] x {0}), one observes that the minimal second eigenvalue 
(as a function of Xq) such that the two nodal domains give rise by symmetry 
to a 3-partition is obtained for xq = 0. 

Doing the same computation on the diagonal, we observe also numerically 
that the minimal second eigenvalue (as a function of the point on the di- 
agonal) such that the two nodal domains give rise by symmetry to a three 
partition is obtained at the center. 

Admitting these two numerical results, what we have proved for the square 
is that the two minimal second eigenvalues are equal. 
This could also suggest that we have a continuum of minimal 3-partitions. 
This point is not completely clear but could result of the analysis of the sin- 
gularity at {0}. 
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When minimizing over Dirichlet-Neumann-Dirichlet or Neumann-Dirichlet- 
Neumann problems in Q'-^, the numerical computations (see [9l §3.3]) suggest 
that the nodal sets of the second eigenfunction never create a 2-partition of 
fl'''^ leading by symmetry to a 3-partition of fl. The corresponding eigen- 
modes lead to too high energy, hence do not qualify as possible candidates 
for minimal 3-partitions. 

Remark 6.1 

Note that we assumed that the minimal partition is symmetric with respect to 
an axis of symmetry; the numerical experiments make this assumption plau- 
sible. One cannot a priori exclude that the first excited K-real eigenfunction 
of the Aharonov-Bohm hamiltonian consists of one line joining to dQ and 
another line joining in Cl two points of dVt. 



6.2 The symmetries of the square 

We now consider a convex domain which in addition to the invariance by (Xi 
and (72 has an invariance by rotation (centered at the origin) of |. We 
have typically in mind the case of the square. This rotation can be quantized 
by 

7^50=^r_.■), (6.1) 



where Tq, is the rotation by a in the plane. We observe from (5.6) that this 
quantization commutes with the operator: 

AA7^| = 7^|AA , (6.2) 

(and with its Dirichlet realization in Cl). 
More generally, we have the following Lemma 

Lemma 6.2 

If u is a K-real eigenfunction of the Aharonov-Bohm Hamiltonian, then u 
and e^^lZiLU are linearly independent K-real eigenf unctions. 

Proof: Let us first verify the K-reality. 

We note that 

TT 

TZk Kv = TZk expiO Tv = exp ~i— exp i6T TZkv , 

2 2 2 ^ 

hence 

TT 

TZe. K = exp —i— KTZiL . 

2 
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This can be rewritten in the form 

TX \ 

exp z— TZiL I o K 



Ko 



TT 



exp z— TZe. 
^ 4 2 



(6.3) 



This proves the first statement. 

We now show that exp i| and u are hnearly independent (over M) inside 
the real space L\ . Let us look at the points of the nodal set belonging to 
the exterior boundary. Their cardinality should be odd by a result of ^IT| 
on K-real eigenf unctions. If u and e^^lZ^u were proportional, this subset 
should be invariant by rotation of | and should have consequently an even 
cardinality. Hence a contradiction. □ 



Proposition 6.3 

In the case of a convex domain having the symmetries of the square, the 
Dirichlet-Neumann problem for the Laplacian on the four half-domains re- 
spectively defined by 



Q+-dh 



Qn{x 
Qn{x 



y <0} , 
y>0} , 



Qn{x 
Qn{x 



y>o}, 
y<o}, 



are also isospectral to the problems introduced in Proposition\5.4 



(6.4) 



K 

Neumann 



n- 





Figure 6: Domains fl 



-dh 



^++dh ^—+dh — dh 



Proof: We explain below how to get the proposition (which has also 
been verified numerically). We start from v G L'](y,2- consider 

.TT 

w = V + exp2— iZnv . 

4. 2 

We have already shown that w is not zero and hence is an eigenfunction. It 
remains to analyze its zero set which should contain an half-diagonal. 
Let us introduce 

S^ = exp^^7^^S^ 
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Using the property that 
we can verify that 

= w , Kw = w . (6.5) 
For = and x + iy = rexpiO, we get 

_ TT 37r_ 
wexpz— = w , exp —i—w = w , 

hence w = for 9 = . So the restriction of w to Q ^'^^ multiphed by a 
phase factor leads to an eigenfunction of the DN-problem for the Laplacian 
in 

The converse does not introduce new problems. □ 



6.3 The covering approach 

As in [i9\, we can also rewrite all the proofs by lifting the problem on the 
double covering fl^, using the correspondence between L|- and the real sub- 
space of the functions in L'^{Q^) such that = — -u where S is associated 
with the map a by (T,u){u!) = u{a{uj)). We recall with the notation of [20] 
that a is defined by associating to each point u of Ci^ the other point (t{uj) 
of Cf^ which has the same projection on Cl. Our initial proof was actually 
written in this way but we prefer to present in this paper another point of 
view. 



In a recent paper, D. Jakobson, M. Levitin, N. Nadirashvili, I. Polterovich 
|24j obtain also nice isospectrality results involving Dirichlet-Neumann prob- 
lems. They actually propose three different proofs of their isospectrality 
results. The second one is a double covering argument which is quite close 
to what we mentioned in the previous paragraph. But there is no magnetic 
version and our magnetic examples seem to be new. So it would be interest- 
ing whether the magnetic approach can produce some isospectrality result 
which differ from the class of results in this paper and the papers by Levitin, 
Parnovski and Polterovich [25j| or [27] . 

One should also mention that for these questions of isospectrality a cover- 
ing argument was already present in earlier works of Berard [3 IH E] , Berard- 
Besson [_6j, Sunada [28].... 
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7 Some heuristics on the deformation of sym- 
metric minimal partitions and numerical 
computations 



7.1 Some heuristics 

Here we discuss very heuristically in which general context the numerical 
computations for the family of rectangles can be done. 

One might investigate the special situation where tight upper and lower 



bounds to £3 are available. We recall from (2.7) that 

If we have a family of domains depending analytically on a parameter a, 
such that for some k, 

limAfc(a) = Lk{0) and Xk{a) < Lk{a) for a > , (7.1) 

aj,0 

then we are led to the question how a minimal /c-partition of of Q{a) 

behaves as a tends to zero. In fact we might investigate more directly 2)^(0) . 
For a = the only 2)^(0) are the one which are nodal partitions associated 
to Afc(O) = Lfc(O) . 

We have investigated this situation for rectangles. One has seen that in 
this case the nodal partition can only have critical points at the boundary. 



Let us start with a slightly more general situation and consider a family 
of simply connected domains which depends analytically upon a param- 
eter (3 > 0. We assume that the spectrum and eigenfunctions corresponding 
to H{Q{f3)) have for small (3 > the following properties 

Ai(/5) < X2i/3) < XsiP) < A4(/3) = , (7.2) 

with 

XsiP) < XiiP) for /? > and A3(0) = A4(0). (7.3) 

Here as usual H{Q{f3)) is just —A with a Dirichlet boundary condition. We 
further assume that for < /5 < Po, niu^^jd)) = 2, but that for /5 = the 
eigenspace U3 of Ai(0) contains an eigenfunction u eU^ with /i(-u) = 3. We 
hence have for /? > that X^{I3) is not Courant-sharp and therefore there is 
a £3(/3) > X^{(3). But for /3 = we have £3(0) = A3(0). 
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According to Proposition |2.5 we have for (3 > three types of non bipar- 



tite partitions, (a), (b), (c). But we have observed that for i23(0), there is 
no 5)3(0) which is not bipartite. So we would hke to understand how a non 
bipartite partition T)3{f]) can be deformed so that it becomes bipartite. We 
emphasize that at the moment we have no mathematical tools permitting to 
rigorously prove the validity of these "deformation argument" but numerical 
tests show that they are rather good for predicting what is observed. 

(a) If the family is of type (a), what seems natural to imagine (see Fig- 
ure [7| is that the critical point should move to one point of the boundary 
and that (at least) two lines will start from this point and end at two other 
points of the boundary. These three points are not necessarily distincts. 



Figure 7: Deformation for type (a). 



(b) If the family is of type (b), then the two critical points should again 
tend to the boundary (see Figure [s]). We note indeed that in the case when 
2:1, 2:2 tend to a point z & Q then we get 4 nodal domains. 



Figure 8: Deformation for type (b). 



(c) In the case (c), a new situtation can occur when the two critical points 
tend to one point in Q. One could indeed imagine a deformation of type (c) 
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minimal partitions on a 3-partition which is diffeomorphic to the figure eight 
(see Figure [9]). 





Figure 9: Deformation for type (c). 



Let us recall that in the case of the rectangle we have explicitly verified 
that a limiting minimal partition cannot have a critical point inside the 
rectangle. 



7.2 Numerics in the case for a family of rectangles 

We look at the case when the parameter (3 introduced in the heuristic sub- 
section is the e of the computations and we go from e = 1^3/8 to e = 1. 
We assume that all the minimal partitions are symmetric with respect to the 
horizontal axis {y = 0} and of type (a). This permits to use the argument 
of reduction to an half-rectangle and to use the approach of the Dirichlet- 



Neumann for each value of e. Figures 10 present the evolution of the candi- 



date to be minimal 3-partition for rectangles Rtt^.tt with e from a/3/8 to 1. 
As mentioned in |9], numerical simulations on the half-square with mixed con- 
dition Dirichlet-Neumann-Dirichlet or Neumann-Dirichlet-Neumann never 
produce any 3-partitions of type (b) or (c). 
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(a) i = (h) t^2 (c) t = 3.45 (d) t = 7 
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